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Recent works have demonstrated that the Floquet-Bloch bands of periodically-driven systems 
feature a richer topological structure than their non-driven counterparts. The additional structure 
in the driven case arises from the periodicity of quasienergy, the energy-like quantity that defines 
the spectrum of a periodically-driven system. Here we develop a new paradigm for the topological 
classification of Floquet-Bloch bands, based on the time-dependent spectrum of the driven system’s 
evolution operator throughout one driving period. Specifically, we show that this spectrum may 
host topologically-protected degeneracies at intermediate times, which control the topology of the 
Floquet bands of the full driving cycle. This approach provides a natural framework for incorporating 
the role of symmetries, enabling a unified and complete classification of Floquet-Bloch bands and 
yielding new insight into the topological features that distinguish driven and non-driven systems. 


After the discoverjff^^ and explanatiorP^ of the quan¬ 
tized Hall effects, topology gained new importance as 
a mechanism for generating extremely robust quan¬ 
tum mechanical phenomena. The realization that the 
Bloch bands of solid state systems could possess non¬ 
trivial topological character istics led to the predictiorPHini 
and exp erime ntal discoverjEE^ of whole new classes of 
materialJEm _ the topological insulators and supercon¬ 
ductors - which host a variety of remarkable and po¬ 
tentially useful phenomena. On a th eoretical level, a 
complete topological classificatioiEEl of such systems 
has been developed, predicting a number of new phases. 
However, finding materials that realize these phases re¬ 
mains a very challenging task, with no known examples 
for many topological classes. 

Motivated by the great successes and open challenges 
in the arena of topological matter, many authors have 
begun to explore the possibilities for reali zing t opologi- 
cal phenomena in driven quantum systemal^^n^ZI Time- 
dependent driving offers the opportunity to control a ma¬ 
terial’s properties in a variety of new ways, potentially 
opening new routes fo r studyi ng topological phen omen a 
in solid stat^EI^ atomicP^fi^Iffl^ optical systemPEH^ 

Intriguingly, driven systems may host an even richer 
array of topological phenomena than their non-driven 
counterparts. To date several examples of topologi¬ 
cal phenomena which ca n only be r ealized in driven 
systems have been foundEHMESIH^ such as the ex¬ 
istence of robust chiral edge states in two dimen¬ 
sional systems whose Floquet bands have trivial Chern 
indices^, and pairs of non-degenerate Majorana end 
modes with protected quasienergy splittings in one¬ 
dimensional system^^. This indicates that periodically 
driven systems feature additional topological structure 
beyond that found in non-driven systems. However, a 
unifying principle for understanding and classifying these 
new phenomena remains lacking. 

In this work we answer the question: under what condi¬ 
tions does the evolution of a driven system become topo¬ 
logically distinct from that of a non-driven system? In 
doing so we develop a powerful and general framework 
that can be used to understand the topology of periodi¬ 



FIG. 1. Phase band representation of the evolution operator 
I/(k,t),Eq. (0. a) For a non-driven system, the phase eigen- 
valnes grow linearly in time, b) Here we show phase bands 
of a periodically driven system which are non-degenerate for 
all k and t. The evolution can be smoothly deformed into 
one obtainable in a non-driven system without closing any 
quasienergy gaps, c) In this case the evolution operator fea¬ 
tures non-removable degeneracies which prevent such a defor¬ 
mation. The evolution is therefore topologically distinct from 
that of any non-driven system, d) Illustration of the phase 
band labeling scheme defined in Sec. [I] The shaded region 
indicates the “phase Brillouin zone.” 


cally driven systems. 

In the analysis of periodically driven systems, the Flo¬ 
quet operator, denoted U(T), plays a central role as the 
stroboscopic evolution operator that propagates the sys¬ 
tem forward in time through each complete driving pe¬ 
riod, T. The spectrum of the Floquet operator, given 
by U{T)\Aln) = plays an analogous role to 

the spectrum of the Hamiltonian in a non-driven system, 
with real-valued energies replaced by periodically-defined 
quasienergies, Sn -\-2 'kN/T = for any integer N. For a 
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system on a lattice, the single particle spectrum forms 
bands, the so-called Floquet bands. Throughout this 
work we focus on systems defined on a lattice, with a 
finite number of bands. While knowledge of the Floquet 
bands is sufficient to understand many aspects of the dy¬ 
namics of a driven system, it was recently shown that the 
topological properties of the evolution are in particular 
not described by U{T) alon^^. A proper description of 
the topology of driven systems must take into account 
the full evolution U{t) for times t throughout the entire 
driving period, 0 <t <T. 

As a means of elucidating the nature of the evolution 
C/(t), we focus on the “phase bands” of the system, i.e., 
the time-dependent spectrum of the system’s evolution 
operator throughout one driving period. As depicted in 
Fig.ig for each time t in the interval 0 < t <T, the eigen¬ 
values of the Bloch evolution operator [/(k, t) 

form bands as a function of the crystal momentum k. 
For illustration we use a “repeated zone” representation 
for the phase bands, though the complete spectrum is 
contained within a single phase Brillouin zone of width 
27r, as indicated by the shaded region in Fig. [^. As a 
function of time, these phase bands form sheets which, 
along with the corresponding eigenvectors, contain full 
information about the evolution of the system. 

Below we examine smooth deformations of the phase 
bands which, keeping the Floquet operator U{T) fixed, 
determine when a given system’s evolution can be 
smoothly deformed into one obtainable in a non-driven 
system. At time t = 0 the evolution is the identity. 
Therefore all phase bands must originate with phases 
(j) which are integer multiples of 27r. For a non-driven 
system with Hamiltonian H, the evolution operator is 
given by U{t) = In this case the phase bands 

diverge from one another linearly in time due to the 
linear phase winding (j) = Et for each eigenstate of 
H with energy E (see Fig. [T^) . For the case of a 
driven system as shown in Fig. Ijjo, the phase bands can 
be straightened through a continuous deformation, such 
that the evolution becomes indistinguishable from one 
generated by a time-independent Hamiltonian. Crucially, 
as we show below, phase bands may be connected via 
topologically-protected degeneracies, or “topological sin¬ 
gularities,” which prevent the evolution from being de¬ 
formed into the canonical form for a non-driven system 
(Fig.[^). These singularities play a central role in defin¬ 
ing the topology of periodically driven systems. 

After establishing the existence of topological singu¬ 
larities in the bulk evolution, we study their ramifica¬ 
tions for the edge mode spectrum of U{T) for a two- 
dimensional system defined in a geometry with edges. 
To this end we identify a complete set of independent 
topological quantities which are (by definition) invari¬ 
ant under any smooth deformation of the bulk evolution 
that preserves the Floquet operator U (k, T). Noting that 
the net number of chiral edge modes in each gap must 
also be invariant under such smooth deformations, we 
derive a bulk-edge correspondence in terms of these in¬ 


variant quantities [Eq. ( |Io| ) below]. We then show how 
the method can be extended to systems in arbitrary di¬ 
mensions, also including the role of symmetries, thus pro¬ 
viding means for a complete topological classification of 
Floquet-Bloch systems. We show that symmetries con¬ 
sidered previously, e.g. in Refs. [22l and l44l which gener¬ 
alize the Altland-Zirnbauer symmetry classes to the case 
of periodically driven systems, can be naturally incorpo¬ 
rated into the phase-band picture. Importantly, we find 
that these symmetries can protect new types of topolog¬ 
ical singularities in the bulk. 

In general we find that, for each bulk gap, the edge 
mode spectrum of a driven system in a given symmetry 
class has the same set of protected features as that of 
a non-driven system in the corresponding class. How¬ 
ever, the global edge mode spectrum and the relation be¬ 
tween edge modes and bulk bands can be quite different. 
Examples of such new or “anomalous” edge phenomena 
include Floquet-Majorana edge modes22E3 quasi¬ 
energy tt/T and chiral edge mode^^in a 2D system with 
topologically trivial bulk Floquet bands. Here we also 
show that periodic driving, for example, allows two-band 
systems with time-reversal symmetry to have helical edge 
modes, while a minimum of four bands is required in the 
non-driven case. 

Interestingly, we find that all the above phenomena 
are closely connected with the appearance of topological 
singularities in the bulk evolution. Due to the additional 
freedom presented by time-dependence, we further specu¬ 
late that there may be other new types of symmetry con¬ 
ditions (beyond those familiar from non-driven systems) 
which can protect new types of topological singularities 
and anomalous edge mode phenomena. 

The remainder of the paper is structured as follows. In 
Sec. [T] we formalize the description of phase bands, and 
characterize the singularities which may prevent them 
from being deformed into a trivial configuration. Then 
in Sec. [IT] we cast the topological characterization of two- 
dimensional systems (without symmetries) in terms of 
the phase bands and their singularities, giving new in¬ 
sight into the winding number invariants found previ¬ 
ously in Ref. 021 In Sec. |HI| we show how additional sym¬ 
metries (e.g., time reversal or particle-hole symmetry) 
can be naturally incorporated into this picture through 
their abilities to protect new types of singularities. Fi¬ 
nally, in Sec. |IV| we summarize our results and discuss the 
outlook for future work. Technical aspects of derivations 
are provided in appendices. 


I. PHASE BANDS OF THE EVOLUTION 
OPERATOR 

We now study the question of when the evolution of 
a periodically-driven Floquet-Bloch system can be con¬ 
sidered topologically distinct from that of a non-driven 
system. In order to do this, we begin by defining the 
phase band picture of Floquet-Bloch evolution. In this 
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section we focus on “bulk” systems with discrete trans¬ 
lation symmetry (with infinite extent or periodic bound¬ 
ary conditions). Here, the crystal momentum k is a good 
quantum number. For now we leave the number of spa¬ 
tial dimensions arbitrary. 

The evolution of a periodically-driven quantum sys¬ 
tem may equivalently be prescribed in terms of either 
a Hamiltonian + T) = H{t), where T is the driv¬ 
ing period, or by the corresponding evolution operator 
U{t) = Te“*-^o , where T denotes time ordering. 

In this paper we primarily work directly with the evolu¬ 
tion operator U{t), which most clearly exposes the topo¬ 
logical features of the evolution. Importantly, although 
the Hamiltonian satisfies H{t + T) = Hit), the evolution 
operator U(t) is generally not periodic in time. 

For bulk systems, crystal momentum k and time t 
parametrize a family of Bloch evolution operators 17(k, t), 
which act within the space of periodic Bloch func¬ 
tions. When the time-dependent Hamiltonian is local 
and bounded, U (k, t) is continuous in crystal momentum 
and time. 

As an important first step in our analysis, we express 
U (k, t) in terms of its spectral decomposition 

N 

C/(k,t) = ^P„(k,t)e-*^"('^’‘), (1) 

n—1 

where P„(k, t) is the projector onto the n-th eigenstate 
of t/(k, t) and is the corresponding eigenvalue. 

Here N is the number of bands in the system. 

We refer to the functions {(^ri(k, t)} as the phase bands 
of the system. In contrast to the quasienergy bands asso¬ 
ciated with a driven system’s Floquet operator C/(k, T), 
the phase bands depend on time, and are continuously 
defined throughout an entire driving cycle, 0 < t < T. At 
time t = T, the phase bands coincide with the system’s 
Floquet bands. An illustration of phase bands for a one¬ 
dimensional system with two bands is shown in Fig. 

To resolve the ambiguity of the labeling of eigenstates 
of U (k, t) we now define a prescription for assigning the 
values of the n indices. We focus on the phase bands 
{(/)„(k, t)}, and work in a repeated zone representation 
where the spectrum is copied and shifted through all in¬ 
teger multiples of 2tt. Recalling that f7(k, 0) = 1, each 
phase eigenvalue must start from an integer multiple of 
2tt. However, as mentioned above, the full spectrum 
of U{k,t) is contained within one “phase 
Brillouin zone.” 

While in principle we could choose ((>„ (k, 0) to be equal 
to any integer multiple of 27r, we choose to work in a fun¬ 
damental phase Brillouin zone in which all phases orig¬ 
inate from zero, i.e., ((>„(k,0) = 0, for n = 1...N (see 
bold curves in Fig. Si)- Next, we demand that each 
4>n{k,t) is a real-valued continuous functiorP^ of both 
k and t. Finally, we impose an ordering condition: if 
(j)n{k,t) > t) for one point in k,t-space, then this 

relation must hold for all k, t. By ordering the indices 



FIG. 2. Graphical depiction of the deformation described in 
Sec. [T] where the time-evolution operator of a periodically 
driven system is deformed into the time-evolution operator 
of a non-driven system. This deformation is always possi¬ 
ble if the time-evolution operator has no degeneracies. After 
flattening, the linear ramp region is expanded to the entire 
interval 0 < t < T, and the bands are straightened into a 
form as in Fig. Si- 



such that n > m implies > (jjm, this prescription de¬ 
fines a unique labeling of the phase bands. 

The ordering condition above is constructed such that 
if two phase bands become degenerate at a particular 
value of k,t, the bands do not “cross” (e.g., the index 
m stays with the lower branch everywhere, while the in¬ 
dex n stays with the upper branch). This arrangement is 
in particular maintained when a band in the fundamen¬ 
tal zone meets a band originating from a different zone, 
see for example band #3 in Fig. Si. Such degeneracies 
between phase bands associated with different branches 
play an essential role in defining the topological charac¬ 
teristics that distinguish driven and non-driven systems. 

We now use the phase band picture to demonstrate 
when it is possible to continuously deform the evolution 
of a Floquet-Bloch system into that of a non-driven sys¬ 
tem, while keeping U{k,T) fixed. Naively, the continuity 
of U (k, t) and of the phase bands might lead one to ex¬ 
pect that the projectors Pn{k,t) are continuous as well. 
If this were true, any continuous deformation of the phase 
bands (j)n(k,t) would preserve the continuity of the evo¬ 
lution operator. It would then always be possible to de¬ 
form the evolution into that of a non-driven system us¬ 
ing a two-step “band-flattening” procedure (see Fig. |^. 
First, for every n = 1... N, deform (()„(k, t) to zero for 
all 0 < t < T — i5t until a small time-interval St before T, 
after which it grows linearly to its final value. If the inter¬ 
val St is small enough, we can assume that the projectors 
are constant there, P„(k,t) = P„(k,T). In the second 
step, let St —>■ T, while keeping the projectors constant 
throughout the linear ramp of the phase. The deformed 
evolution is now identical to that of a non-driven system 
with the Hamiltonian 

Mk) = ^5]</-„(k,r)P„(k,T). (2) 

n 

The picture above seems to imply that all periodically- 
driven systems are topologically equivalent to non-driven 
systems (i.e., they can be related by smooth deformations 
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that keep the Floquet operator fixed). However, th e ex- 
istence of phenomena such as anomalous edge modeJ^^EHl 
shows that this cannot be the case. 

Where could the argument break down? In the first 
step, we assumed that the phase bands could be contin¬ 
uously deformed to zero throughout the entire driving 
period, up to a short interval 5t in which the projec¬ 
tors were assumed to be constant. However, in principle 
one may imagine that the evolution operator could host 
degeneracies around which the projectors are discontin¬ 
uous (the degeneracy of the eigenvalues ensures that U 
stays continuous). In the presence of such a discontinu¬ 
ity, the degeneracy could not be lifted without breaking 
the continuity oi U. In this way a phase band in the 
fundamental zone may become “glued” to another band 
from a neighboring branch of phases (see Fig. &)■ 

In the absence of the discontinuities described above, 
the evolution of any driven system can be smoothly de¬ 
formed to that of a non-driven system, as in Fig. and 
anomalous edge states would be impossible. Thus we 
are led to the unavoidable conclusion that the evolution 
operators of periodically driven systems must support 
topologically-protected degeneracies. In the next subsec¬ 
tion we show explicitly that such degeneracies can exist 
in two-dimensional (2D) systems. In Sec. HI we gener¬ 
alize to other dimensions and symmetry classes. Below 
we refer to these topologically-protected degeneracies as 
“topological singularities.” 


A. Topological singularities in two-dimensional 
systems 

In this subsection we explicitly demonstrate the exis¬ 
tence and nature of topological singularities in the evolu¬ 
tion operators of two-dimensional systems. We further¬ 
more show that in a region in k, t-space where U (k, t) 
is degenerate, the degeneracy can either be lifted every¬ 
where or reduced to a cluster of isolated singularities. 

Let U (k, t) be the bulk time evolution operator of a 
two-dimensional system with no symmetries other than 
the discrete translational symmetry of the lattice. Con¬ 
sider now a point Sg = (kg, tg) in k, t-space where two ad¬ 
jacent phase bands, m and m', are degenerate (mod 27r). 
At Sg, the degenerate subspace is spanned by the Bloch 
states liJm) and Due to the continuity of U(k,t) 

and the existence of gaps to other phase bands, we can as¬ 
sume that the subspace spanned by the two intersecting 
bands is constant within some finite sized neighborhood 
around Sg in k, t-space. The remaining non-degenerate 
bands {|Xn)} and their associated phases {</>«} can also 
be assumed to be constant within this neighbourhood. 
Close to Sg, the time evolution operator thus takes the 
form 

U{s)=Y^ IXn)(X«|e"*'^" + 

( 3 ) 


where M is a 2 x 2 unitary matrix and we parametrize 
the three-dimensional (k, t) space by a single variable s. 

The unitarity of M means that we can write it as 

M(s) = exp [-i(t)d{s) - ifj{s)crj] , (4) 

where summation over repeated indices is used. Here 
is a real-valued function, whose value at Sg gives 
the common eigenvalue of the two degenerate bands, {aj} 
are the Pauli matrices, and {/j(s)} are real continuous 
functions that satisfy /i(sg) = 0. 

We assume that U, and thereby /, is differentiable in 
a neighborhood around Sg, and expand fj to linear order 
in (s — Sg) around Sg. Noting that /j(sg) = 0, we write 

M(s) « exp [-i(^d(s) - i(s - So)jSjk(Jk], (5) 

where Sjk = djfk{so) is a real 3x3 matrix. The case 
where the linear term in (s — Sg) also vanishes will be 
covered shortly. 

We first consider the case where the matrix S has 
rank three, such that the coefficients of all three Pauli 
matrices vary independently as s explores the neighbor¬ 
hood around Sg. In this case, the degeneracy is topologi- 
cally protected, similar to the case of a Weyl node^ an 
infinitesimal change of the time evolution operator can 
never lift the degeneracy, but rather can only infinitesi¬ 
mally shift the location where it appears. A single such 
degeneracy can thus not be lifted with a continuous de¬ 
formation of the system, and is therefore topologically 
protected. We thus define a topological singularity of a 
two-dimensional system to be an isolated degeneracy of 
the time evolution operator where the matrix S describ¬ 
ing the linearization of U in its neighborhood [Eq. ([^] is 
invertible. 

In addition to the isolated singularities described 
above, we may also find cases where S is not invertible. 
This occurs when two phase bands are degenerate along 
a line, surface, or three-dimensional region in k, t-space, 
such that Sg is one point on this manifold. In such cases, 
the rank of S is equal to 3 — D, where D is the dimension 
of the degenerate manifold. These extended degenera¬ 
cies are not topologically protected: the degeneracy can 
generically be lifted in a neighborhood of Sg with a local 
perturbation, letting /^(s) —>• /^(s) -|- 6g{s)vi in Eq. 

Here S controls the strength of the perturbation, v is a 
3-dimensional vector satisfying ViSij = 0, and ^(s) is a 
real continuous function that vanishes outside a neigh¬ 
bourhood around Sg, within which can be taken 

to be constant. 

Importantly, the local perturbations described above 
only lift the degeneracy patch-wise, in one small region 
at a time. If one tries to lift the degeneracy over the 
entire manifold, two cases are possible: either the degen¬ 
eracy can be lifted everywhere, or there will be a dis¬ 
crete set of points where topological singularities remain. 
Hence we conclude that, if the time evolution operator is 
degenerate throughout a finite-dimensional manifold, it 
is always possible to apply an infinitesimal perturbation 
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that either completely lifts the degeneracy, or reduces it 
to a cluster of topological singularities. 

With the existence of topological singularities estab¬ 
lished, we now further characterize their properties. Each 
singularity can be assigned a charge (or vorticity) q: 

q = sgn[det5'], (6) 

where S is the linearization of /fe(s) around Sq, see 
Eqs. 0 and ([5|. 

In two-dimensional systems, the charges of topologi¬ 
cal singularities have direct connections with the Chern 
numbers of the phase bands. Consider the “instan¬ 
taneous” Chern number of phase band n, Cn{t) = 

/ (fkTT:{Pn[dk^Pn,dk^Pn]}. As long as no singular¬ 
ities are encountered, the Chern number Cn{t) cannot 
change due to the continuity of Pn(t). However, when 
two phase bands meet at a singularity with vorticity q, 
the Chern number for the “upper” phase band changes 
by q as the singularity is traversed in time, while the 
Chern number of the other band changes by —q. Here, 
the “upper” band is band to -I- 1 if the singularity con¬ 
nects bands to and m -|- 1, and band 1, if the singular¬ 
ity connects band 1 and N through the phase Brillouin 
zone edge. As a consequence of the argument above, any 
driving protocol that yields Floquet bands with differ¬ 
ent Chern numbers from those of the initial Hamiltonian 
H{0) must induce one or more topological singularities 
in U{k,t). 

In this section we showed that the evolution operator 
of a periodically-driven system may host topologically- 
protected degeneracies, or “topological singularities.” As 
we concluded in the beginning of the section, topologi¬ 
cal singularities can obstruct the smooth deformation of 
the evolution of the driven system into that of a non- 
driven system. Specifically, in the case where the “bot¬ 
tom” and “top” phase bands, 1 and N, are connected by 
a singularity through the phase Brillouin zone-edge, their 
respective phase values at the singularity must differ by 
2tt. In this situation it is impossible to simultaneously 
flatten both of the bands to zero (compare Figs. 0, d 
with Fig. 0. In contrast, for singularities that do not 
pass through the phase zone edge (i.e., those connecting 
bands to and to -|- 1, with to < N), the two correspond¬ 
ing phase eigenvalues coincide at the singularity. In this 
case nothing prevents deforming the two phase bands si¬ 
multaneously to zero, thereby removing the singularity. 
Hence singularities of the first type, i.e., “zone-edge sin¬ 
gularities,” are special: it is precisely these singularities 
that cannot be eliminated by smooth phase-band defor¬ 
mations, thus distinguishing driven from non-driven evo¬ 
lution. In Sec. 10 below we formulate the topological 
classification for two-dimensional periodically driven sys¬ 
tems in terms of the phase bands and zone-edge topolog¬ 
ical singularities, and derive the corresponding bulk-edge 
correspondence. 


B. Natural quasienergy zone 

Before moving on to classification, we briefly introduce 
some further labeling notation that will be useful for re¬ 
ferring to specific Floquet bands and gaps in the dis¬ 
cussion below. In particular, we apply the phase band 
labeling prescribed in the beginning of this section to the 
Floquet bands themselves. We label Floquet band n such 
that e„(k) = ^„(k, T)/T (see Fig. 0). This uniquely 
determines “natural” quasienergy band indices for a sys¬ 
tem, and the bands {Sn} for n = 1... A^ define a natural 
quasienergy (or Floquet) Brillouin zone. 

The Floquet band indices are topologically invariant: 
within the predefined labeling scheme, the index of a Flo¬ 
quet band cannot change under any smooth deformation 
that keeps the gaps above and below the band open. No¬ 
tably, the natural indices of the Floquet bands can be 
found directly with the help of the time-averaged Hamil¬ 
tonian, without having to track the full time development 
of the phase bands. In particular, note that 

en{k) = ^£dtTr[H{k,t)]. (7) 

Recall that all en(k) lie within an interval of width 2Tr/T, 
and that the bands are ordered by increasing n. From 
this information and Eq. 0 above, the indexing of Flo¬ 
quet bands is uniquely determined. 

Below we will also apply the labeling scheme to the 
quasienergy gaps. For the following discussion, we refer 
to the quasienergy gap above band to as gap to. Due 
to the periodicity of quasienergy, a driven system with 
N bands has an additional gap (as compared to a non- 
driven system), which separates band N from band 1, 
across the quasienergy zone edge. We thus refer to gap 
N as the “zone-edge gap” of the Floquet spectrum. In 
the non-driven limit T —0, the zone-edge gap becomes 
infinitely wide, while the other gaps remain finite. 

II. TOPOLOGICAL CLASSIFICATION OF 
FLOQUET-BLOCH SYSTEMS IN TWO 
DIMENSIONS 

Having introduced the concept of phase bands and 
demonstrated the existence of topological singularities, 
we now consider the implications of these results for 
the topological properties of a two-dimensional Floquet- 
Bloch system with no symmetries. By “topological” we 
mean those properties that are invariant under any con¬ 
tinuous deformation of the bulk time-evolution operator 
that preserves its continuity in crystal momentum and 
time, and keeps the bulk gaps open in the quasienergy 
spectrum of the Floquet operator U{T). Any such quan¬ 
tity is a topological invariant of the system. Importantly, 
this definition means that topological invariants must be 
independent of the choice of time origirP^. 

Analogous to a Chern insulator, a two-dimensional 
Floquet-Bloch system defined in a geometry with edges 
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FIG. 3. Continuous deformation of a generic time evolution 
operator which only preserves information about the topo¬ 
logical singularities spanning the zone-edge gap (region I), 
and the Floquet bands (region II). These features cannot be 
removed, since the Chern indices of the Floquet bands and 
the total charge of the zone-edge singularities can only be 
changed by closing at least one quasienergy gap. Any topo¬ 
logical invariant should be expressible in terms of only these 
non-removable features. Note that the singularity in the gap 
around phase (f> = 0 is eliminated under the deformation. 


may host protected chiral edge modes within its bulk 
quasienergy gaps. The chiral edge modes are topologi¬ 
cally protected, meaning that the net number nedge(^7i) 
of chiral edge modes in bulk quasienergy gap m is invari¬ 
ant under continuous deformations of the bulk evolution 
operator 17 (k, t), or equivalently of the Hamiltonian H{t), 
that keep quasienergy gap m open. Thus we recognize 
tiedge(w) as a topological invariant of the system. 

In this section we demonstrate that considerations 
about the bulk phase bands allow us to identify all in¬ 
dependent topological invariants of a two-dimensional 
Floquet-Bloch system. Subsequently, we use these invari¬ 
ants to construct the bulk-edge correspondence, provid¬ 
ing a direct link between the edge mode spectrum and the 
bulk phase band properties of two-dimensional Floquet- 
Bloch systems. 


A. Topological invariants of 2D systems 

In this subsection we identify the properties of a two- 
dimensional Floquet-Bloch system that are invariant un¬ 
der smooth deformations of the phase bands. In Sec. 
we found that the phase bands of a 2D system can safely 
be deformed anywhere in k, t-space except at isolated 
topological singularities. Therefore we may expect that 
the singularities play an important role in the topological 
classification of periodically driven systems. 

To elucidate the importance of topological singularities 
we consider the following deformation of a time-evolution 
operator 17(k, t), shown schematically in Fig. With¬ 
out changing the eigenstate projectors P„(k,tj7 deform 
the phase bands (/)„(k,t) to zero everywhere except for 
within small isolated regions that surround each zone- 
edge singularity and within a small time-interval St be¬ 


fore T, where the phase bands wind linearly to their fi¬ 
nal values. These final values, which define the Floquet 
bands, are kept fixed under the deformation. We refer 
to the region t < T — St in k,t space as region I, and 
to the final region t > T — St as region II (see Fig. |^. 
Any time-evolution operator can be deformed in this way 
such that continuity is preserved and no quasienergy gap 
is closed. Without changing any topological invariant, 
the deformation effectively discards all information about 
the time-evolution operator except for the phase bands 
at time t = T (i.e., the Floquet bands themselves), and 
the zone-edge singularities. 

Consider now the remaining features of the phase band 
structure that could not be smoothly deformed away. We 
found in Sec. m that it was possible to change the loca¬ 
tion k, t of each singularity through a continuous defor¬ 
mation. Through such a deformation, it is furthermore 
possible to create and annihilate pairs of zone-edge singu¬ 
larities with opposite charges. Hence the only invariant 
quantity we can associate with region I is the sum of the 
charges of all zone-edge singularities, J2i 

For region H, we note that at t = T any two projectors 
P(k,T) and P'(k, T) can be continuously deformed into 
each other if and only if their Chern numbers are the 
sam^. Hence the only independent invariants we can 
associate with region H are the Chern numbers of the 
individual phase bands at t = T. 

The arguments above show that a two-dimensional 
Floquet-Bloch system with N bands has exactly N inde¬ 
pendent topological invariants characterizing it. These 
invariants are the integers 



where C„ is the Chern number of Floquet band n (see 
Sec. IB for definition of the quasienergy band indices). 
The index i in the sum runs over all topological singu¬ 
larities in the zone-edge gap. The Chern number of the 
last band Cn is not included since Cn = 0. We 
see that while an fV-band non-driven system is charac¬ 
terized by — 1 independent integer-valued (Z) invari¬ 
ants (the Chern numbers of each of the N — 1 lowest 
bands), Floquet-Bloch systems are characterized by N 
integer (Z) topological invariants. The additional invari¬ 
ant is the net charge of the topological singularities in 
the zone-edge gap. 


1. Bulk-edge correspondence for two-dimensional 
Floquet-Bloch systems 

We now seek to derive a bulk-edge correspondence that 
gives the net number of chiral edge states that will appear 
within a given gap m of the bulk Floquet spectrum when 
the system is defined in a geometry with an edge. To 
this end we seek the combinations of the N numbers in 
Eq. Q that remain invariant when all gaps except for 
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gap m are allowed to close (see Sec. IB for the labeling 
convention for the quasienergy gaps). We know at least 
one such quantity should exist since the number rie^ge{m) 
of edge modes in gap m has this property. 

In order to find the invariant combinations, we note 
that the Chern numbers of the individual bands I to m 
can be changed by closing the quasienergy gaps between 
them. Only their sum Sm = remains constant 

under such operations. Furthermore, if m 7 ^ IV, all zone- 
edge singularities can be removed through the plane t = 
T by closing the zone-edge quasienergy gap (i.e., the gap 
between band N and band I). Importantly, however, 
the Chern number of band I, and thereby Sm, changes 
by q each time a singularity of charge q is removed in 
this way (see the discussion at the end of Sec. IA on the 
relationship between Chern numbers and singularities). 
Hence there only exists one combination of the invariants 
in Eq. ([^ which remains invariant when only gap m is 
required to stay open: 


( 9 ) 

n—1 i 

Any two evolutions characterized by the same value of 
the invariant Wm can be smoothly deformed into one an¬ 
other without closing quasienergy gap m. Crucially, this 
tells us that the number of chiral edge modes riedgeinT,) 
in gap m should be some function of Wm, and possi¬ 
bly m itself. Standard spectral flow arguments require 
^edge(w) —'nedge(w— 1) = Cm- The Only way this can be 
realized is if nedge(™) = Wm + K for some universal con¬ 
stant K. Considering the trivial special case H{t) — 0, 
where both wn and UedgeiN) are zero, we find that K 
must be zero. We thus arrive at the following new re¬ 
sult for the net number of chiral edge modes in a two- 
dimensional system: 


riedgeim) ■ (10) 

n—1 i 


The simple expression above provides a direct way of 
evaluating the edge mode count given by the winding 
number formula found in Ref. 14111 The first term is the 
result one obtains simply when analyzing a non-driven 
system with the phase band framework, taking T to be so 
small that the phase bands do not cross. The second term 
has no equivalent in non-driven systems, and accounts for 
the anomalous edge modes that were discussed in Refs. HU 
and 031 Additionally, Eq. (10 1 shows that the number 
of edge modes in the zone-edge gap is given by the net 
charge of all zone-edge topological singularities. 


In Appendi x[A| we provide an explicit derivation show¬ 
ing that Eq. ( |10[ ) is equivalent to the winding number 
formula of Ref. 1431 Below we refer to Wm\U] as the wind¬ 
ing number of U in gap m. 


B. Topological singularities in a specific 2-band 
model 


To make our discussion more concrete, in this sub¬ 
section we demonstrate the results above on a variation 
of the explicit model considered in Ref. 031 Consider 
a tight-binding model on a 2D bipartite square lattice, 
described by the time-dependent Bloch Hamiltonian 


i7(k,t) = V Jnit) (cr+e' 






)+Va,, ( 11 ) 


where az and cr^ = (cr^ ± iay)/2 are the Pauli matrices 
acting in the sublattice space, and the vectors {b„} are 
given by bi = —b 3 = (a, 0 ), and b 2 = —b 4 = ( 0 , a), with 
a being the lattice constant. In real-space, Hamiltonian 
0 consists of hopping terms between nearest neigh¬ 
bour sites on the bipartite lattice. The Hamiltonian is 
T-periodic in time. Each driving cycle consists of five 
time intervals of length T/5, with J„(f) = Xn during the 
n-th interval, while all the other hopping amplitudes are 
set to zero. In the fifth interval, all hopping amplitudes 
are zero while the sublattice potential V remains on. 

In Ref. 031 anomalous edge modes were observed in 
the case where A„ = J, for certain ranges of the pa¬ 
rameters J and V. According to the discussion in the 
subsections above, this implies that topological singular¬ 
ities are present. Indeed, when in a nontrivial phase, the 
two phase bands touch through the zone-edge along the 
line kx = ky, at a specific time that depends on param¬ 
eter values. To demonstrate that this degenerate region 
contains topological singularities, we add a small time- 
dependent perturbation to break the extended deg ener- 
acy into isolated singular points (see section I A). We 


implement the perturbation by reducing the hopping in 
the y-direction slightly compared to the cc-direction, such 
that Ai = A3 = J and A2 = A4 = (I — a)J, where a is 
a small parameter. We then numerically calculate the 
time-evolution operator at a representative set of points 
in (k, t)-space for the parameter choice J = —2.5tt/T, 
V = O.Stt/T, and a = 0.2. From diagonalization of the 
time-evolution operator we obtain the phase band struc¬ 
ture of the model, and hnd four topological singularities 
(see Fig.|^. One singularity has charge —1 and connects 
the two bands through the zone edge, while the other 
three have charges 1 ,- 1 , and — 1 , but do not cross the 
zone-edge. The charges are found numerically. 

In Figs. |4^-c, the phase band structure is plotted for 
three values of fixed k^. The k^ values are chosen where 
the four topological singularities appear (two of the sin¬ 
gularities appear at the same kx). The Chern numbers 
of the Floquet bands are zero. 

Next we confine the model to a strip geometry with 
edges parallel to the y-direction, by truncating the 
real-space Hamiltonian of the model in the ^-direction. 
We numerically calculate the Floquet operator of this 
truncated tight-binding Hamiltonian and obtain the 
quasienergy band structure shown in Fig. 1^. On each 







FIG. 4. Explicit demonstration of topological singularities and anomalous edge states, a-c) Phase band structures for the 
model in Eq. (Ill for fixed values of k^, with singularity charges (all ±1) indicated. The values are a) 0.245/a, b) 1.533/a, 
c) 2.084/a. d) Quasienergy band structure of the model in a strip geometry. Both bands have Chern number zero, and we find 
one chiral mode on each edge, in each quasienergy gap. Edge modes on opposite edges are indicated by different colors. 


edge we find the net number of chiral edge modes to be 
1, in both bulk quasienergy gaps. This behavior is fully 
consistent with result (10) above. 


III. TOPOLOGICAL CLASSIFICATION OF 
FLOQUET-BLOCH SYSTEMS WITH 
SYMMETRIES 

In the previous section we showed that the richer topo¬ 
logical structure of two-dimensional periodically driven 
systems was due to the possibility of non-removable sin¬ 
gularities in the phase bands of such systems. Building 
on this result, we now seek to describe how additional 
restrictions on the evolution (e.g., as imposed by discrete 
symmetries) can protect new types of phase band singu¬ 
larities in one, two, or three dimensions. These new sin¬ 
gularities provide the basis for a symmetry-based topo¬ 
logical classification of Floquet-Bloch systems. 

Inspired by the rich structure of the perio dic t able 
of topological insulators in non-driven systemj- ^^ l ^^ ( we 
focus on driven system analogues of the ten Altland- 
Zirnbauer (AZ) symmetry classes. In the first subsection 
below we describe two types of symmetry conditions on 
the evolution operator (“instantaneous” or “time non¬ 
local”) which provide useful ways of generalizing the AZ 
symmetries to driven systems. 

Note that the instantaneous and time non-local condi¬ 
tions are chosen as illustrative examples to demonstrate 
the power and adaptability of the phase band framework. 
These conditions are not necessarily the only ways of gen¬ 
eralizing the AZ symmetries. Note also that these are not 
the only types of conditions that can protect singulari¬ 
ties - it will be an interesting direction for future work to 
seek wholly new types of symmetries on U (t) which may 
protect additional types of topological singularities. 

Following the discussion of smooth phase band defor¬ 
mations from Sec.|^ we find the exhaustive classification 
for one-dimensional systems with particle-hole symmetry 
and identify the related bulk-edge correspondences. We 
then go on to find the bulk-edge correspondences for two- 
and three-dimensional systems with time-reversal sym¬ 
metry. The bulk-edge correspondences that we obtain in 


this section for one- and two-dimensional systems coin¬ 
cide with those found in Refs. and SH respectively. 
To our knowledge, analogous results for the three dimen¬ 
sional case have not been derived before. 

Within the symmetry framework that we consider, we 
find that the edge mode spectrum in a single bulk Floquet 
gap has the same classification for driven and non-driven 
systems. However, the global edge mode spectrum of a 
periodically driven system generally has a richer clas¬ 
sification than its non-driven counterpart, with a cor¬ 
respondingly richer mathematical relationship between 
bulk and edge properties. For each symmetry class, we 
thus find a larger number of distinct topological phases 
than in the corresponding non-driven cases. Interest¬ 
ingly, we find for example that a periodically driven two- 
hand system with time-reversal symmetry can host topo¬ 
logically protected edge modes, while a minimum of four 
bands is necessary in the non-driven case. This richer 
variety of topological phases in driven systems originates 
from the periodicity of quasienergy, i.e., the presence of 
the zone-edge gap in the Floquet spectrum, and the ex¬ 
istence of symmetry-protected topological singularities 
which may reside in the corresponding phase band gap. 


A. Symmetries in periodically driven systems 

We now identify two symmetry conditions on the evo¬ 
lution operator which can protect new types of singular¬ 
ities in the “time-bulk” of the phase band structure (i.e., 
singularities occurring for intermediate times 0 < t < T): 
“instantaneous” symmetries of the form 

u{t) = SU{t)S-^, ( 12 ) 

and “time non-local” symmetries of the form 

U{t) = SU{t* (13) 

Here S may be a unitary or an anti-unitary operator, and 
denotes a special point in the driving cycle. In the sub¬ 
sequent discussion, we always pick the time origin such 
that = T. Note that the instantaneous symmetries 
relate the evolution operator to itself at a given time, 
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while the time-non-local symmetries relate the evolution 
operator to itself at different times. 

The conditions in Eqs. (12) and (13) can be used 
to ensure that the Floquet operator t7(T) and/or the 
corresponding “effective Hamiltonian” HeS, defined via 
U{T) = falls into any one of the ten Altland- 

Zirnbauer symmetry classes. For example, particle-hole 
symmetry is guaranteed via an instantaneous condition 
of the form (12), with S anti-unitary (see below). On 
the other hand, time-reversal and chiral symmetries are 
imposed via time-non-local conditions of the form (13), 
with S anti-unitary and unitary, respectively. Further¬ 
more, just as in the non-driven case, the symmetry con¬ 
ditions where S is anti-unitary divide into two subclasses, 
depending on whether S squares to 1 or —1. 


B. Particle-hole symmetry with = 1 


In this subsection we use the phase band framework to 
develop a topological classification for periodically driven 
systems with particle-hole symmetry (PHS). Here we im¬ 
pose PHS via an instantaneous symmetry condition as 
in Eq. (12), where the operator S is anti-unitary and 
squares to 1 (analogous to symmetry class D in the AZ 
conventiorP^ . Such a condition is naturally satisfied, for 
example, in the Bogoliubov-de Gennes Hamiltonian of 
a driven spinless superconductor. The one-dimensional 
(ID) case was considered previously in Ref. [22l Here 
we use the phase band framework to identify topological 
invariants and to derive the bulk-edge correspondence, 
obtaining results consistent with the findings of Ref. [521 
The condition above implies that “class D type PHS” 
is present if and only if there exists a basis where the 
instantaneous Bloch Hamiltonian of the driven system, 
/i(k, t), satisfies /i(k, f) = —h{—'k,t)*. Consequently, in 
this basis, the evolution operator at each time t satis¬ 
fies [/(k,t) = C/*(—k,t). This furthermore implies that 
the time evolution operator of a particle-hole symmetric 
system with 2N bands can be written in the form 


N 

U{k,t) = [Pn(k,t)e-*^"(‘^’*) +P^(k,t)e-*^-('^’‘)] , 

n=l 

(14) 

where the phase bands {(/„(k, t), 4>n(k., t)} are continuous 
and non-crossing (as defined in Sec.|l]), and, in the basis 
specified above. 


Pn{k,t) = P*{-'k,t), =-(j)n{-k,t). (15) 

The ambiguity of the labeling of bands is removed by 
requiring 0 < (/„ < tt, with the conjugate phases sat¬ 
isfying —TT < (j)n < 0- In previous sections we labeled 
the bands according to increasing n, starting from the 
lowest band. Here, making use of the symmetry of the 
spectrum, we start the labeling from the band with the 
smallest positive phase. 


Similar to the approach in Secs. [l| and [Til we now seek 
to identify topological invariants by considering quanti¬ 
ties that do not change under smooth deformations of 
the evolution which preserve the particle-hole symmetry^ 
expressed via Eqs. (14) and (15). Preservation of the 
symmetry can be ensured by continuously deforming half 
of the phase bands and projectors, {(/„(k, t), P„(k,t)}, 
with their conjugate partners following in accordance 
with condition (15). However, as we found for the two- 
dimensional case considered in the previous section, free 
deformation of the phase bands may be obstructed at 
certain isolated points in k, t-space where protected sin¬ 
gularities are encountered. 

In order to see how a singularity may be protected 
by particle-hole symmetry (in any dimension), consider 
an inversion invariant point kinv in the Brillouin zone, 
where Eq. (15) directly relates the phases and projectors 
of conjugate partner bands. Suppose that at some time to 
two conjugate phase bands m and fh become degenerate 
at kinv) with the labeling we prescribed earlier, m may 
be either 1 or N. Because these bands form a conjugate 
pair, the degeneracy may only occur at phase 0 or tt, see 
Eq. (15). Eor times t close to to, the evolution operator 
at kinv can be written 


[/(ki„v,t)= ^ (lx«)(Xn|e-*^"(‘) + |x:)(x:|e“^'*(‘)) 

n^m 

+ ^ Mab{t)\lpa){i^b\, (16) 

a,6=1,2 


where {|Xn), IXn)} the eigenstates of C/(kinv,to) that 
do not become degenerate (assuming the system has more 
than two bands), and { 11 / 1 , 2 )} are two states that to¬ 
gether span the subspace of the degenerate eigenstates 
m and m at t = to- The 2x2 matrix Jlf(t) is unitary 
and depends continuously on t. Furthermore it satisfies 
Mab{to) = ±/ab, where the sign depends on whether the 
bands meet at phase tt (for —) or at phase 0 (giving -k). 

The symmetry (HI implies that {7(kinv,t) must be 
real for all t. Furthermore, we can take \ip 2 ) to be 
real, since the two degenerate eigenstates {|Xm),IXrn)} 
are complex conjugates of each other. The first sum in 
Eq. (16) is also real, thus implying that M{t) must itself 
be real. As a result, we can write M{t) as 


M{f) = ±e 


- _j_ Iq) 


(17) 


where the parameter A is real. 

The expression for M in Eq. 0 directly manifests the 
fact that the degeneracy is topologically protected: any 
local smooth deformation of the time-evolution operator 
can only continuously change M (t) via the parameters A 
and to, neither of which lifts the degeneracy. The two 
possible signs of M{to) indicate that there can in general 
be two types of singularities at each inversion invariant 
point, namely, singularities occurring at phase 0 and at 
phase TT. 
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Classification of one-dimensional systems with PHS 

We now use the PHS-protected topological singulari¬ 
ties described above to construct the topological classifi¬ 
cation and bulk-edge correspondence for one-dimensional 
systems with particle-hole symmetry. In terms of the 
Floquet spectrum, a system with PHS may host topolog¬ 
ically protected edge modes in its bulk gaps at quasiener¬ 
gies 0 and tt/T. The parities j^ph(O) and i^PH(7r/T) of the 
numbers of edge modes in these two Floquet gaps are 
topological invariants. In P ef. [22l i^ph(O) and j^PH(7r/T) 
were identified as the numbers of times the phase-bands 
at inversion-invariant momenta cross phase 0 and tt, re¬ 
spectively. We now use the framework developed above 
to develop a more complete understanding of this rela¬ 
tionship and to explicitly prove it. 

As for the two-dimensional systems discussed in 
Sec.[Il^, the phase bands of the ID particle-hole symmet¬ 
ric system can be freely deformed (while maintaining the 
symmetry as described above) anywhere except at topo¬ 
logical singularities. Analogous to the procedure depicted 
in Fig. 1^ we deform the phases (j)n for n = 1,..., A^ to 
zero everywhere, except around the phase-7r-singularities 
and in a short interval at the end of the driving pe¬ 
riod where we let the phases wind to their final values. 
This is done in a way that keeps U{k, T) hxed, and such 
that the conjugate bands follow the deformation to pre¬ 
serve particle-hole symmetry via relation (15). Thus it 


is evident that the topological classification of the one¬ 
dimensional particle-hole symmetric system should de¬ 
pend only on the properties of the evolution operator 
around any phase-7r (i.e., zone edge) singularities, and at 
t = T (i.e., on the Floquet bands themselves). 

Which characteristics of the singularities are topolog¬ 
ically protected? For one-dimensional systems, all (zone 
edge) singularities are topologically identical if the sys¬ 
tem has more than two bands: the evolution can be de¬ 


formed such that the vectors IV’ 1 , 2 ) in Eq. (16) are the 
same for all singularities. To see this, note that any two 
pairs of real orthogonal vectors in a complex space of 
more than two dimensions can be continuously rotated 
into each other, if the space has dimension 3 or more. 

When the vectors |'(/'i) and \'ip 2 ) are the same for all 
singularities, it is possible to create or annihilate pairs 
of singularities with opposite sign of the parameter A 
through continuous deformations of the evolution opera¬ 
tor. However, since all (zone edge) singularities are topo¬ 
logically identical as argued above, it is possible through 
a local smooth deformation to interchange {tpi) and \tjj 2 ) 
for an individual singularity, thereby flipping the sign of 
A. Hence all zone edge singularities are identical, and can 
be created or annihilated pairwise through smooth phase- 
band deformations. Thus only the parities Qo and Q.,r/a 
of the numbers of zone edge singularities at fc = 0 and 
fc = 7r/a, respectively, are topologically invariant (here a 
is the lattice constant). The special case of two bands is 
discussed at the end of the section. 

In addition to describing the singularities, we must also 


keep track of the invariants of the Floquet bands, which 
are encoded in the effective Hamiltonian 

Efeff(fc) = ;|,^[(^„(A:,T)P„(A:,T)+<^^(fc,T)P^(fc,T)]. 


The Floquet bands are completely characterized by the 
two standard invariants 770 and 77 ,^ for a non-driven 
particle-hole symmetric system (see, e.g., Ref. HID, where 
r]k = sgnPf [i?eff(fc)] for fc = 0,7r/a. 

Collecting the invariants identihed above, we thus find 
that a generic one-dimensional (translationally-invariant) 
driven system with particle-hole symmetry is fully char¬ 
acterized by the four Z 2 invariants 




(18) 


Compared to the case without driving, translationally 
invariant periodically-driven systems with particle-hole 
symmetry are characterized by two additional Z 2 invari¬ 
ants. The additional invariants relate to phase band sin¬ 
gularities in the zone-edge gap (i.e., at phase tt). 

We now seek a bulk-edge correspondence which gives 
the edge mode parities t'pH(O) and PpH(7r/T) for the gaps 
at quasienergies 0 and tt/T in terms of the four num¬ 
bers in Eq. (18). To see which bulk invariants deter¬ 


mine the edge mode parity ppH( 7 r/T) in the quasienergy 
gap at tt/T, we first note that the bulk evolution can be 
smoothly deformed to “flatten” the Floquet bands to zero 
quasienergy. That is, we may continuously transform 
U{T) —)• 1 (or, equivalently, Hgff —>■ 0), without clos¬ 
ing the gap at tt/T. Such a deformation cannot change 
the parity of the number of edge modes in the open gap, 
and therefore r'pH( 7 r/r) should not depend on 779 or 77,r/a- 
Furthermore, note that a real-space gauge transformation 
can be performed on the ID system to make fc —>■ k-\-TT /a. 
However, a gauge transformation cannot change the num¬ 
ber of edge modes t'pH( 7 r/r) which appear in the gap 
when an edge is created. Therefore Qo ^tnd Q^^/a should 
appear symmetrically in the final expression. Thus we 
are led to the expression: 7 /pH( 7 r/r) = ±(—IjQo-i-QTr/a^ 
Considering the special case H{k,t) = 0 fixes the sign, 
giving 


i^pk{7^/T) = (-l)Qo+Qv», 


(19) 


where ppH(7r/T) = —1 corresponds to an odd number of 
Floquet-Majorana edge modes with quasienergy tt/T. 

For the gap at quasienergy zero, we identify which com¬ 
binations of the numbers in Eq. (18) are invariant when 
we allow the gap at quasienergy ttJT to close. By closing 
this gap, we can remove the zone-edge singularities with¬ 
out changing the number of edge modes at quasienergy 0. 
However, each time we remove a singularity at crystal- 
momentum fc in this way, rjk changes its sign (see Ap¬ 
pendix . Hence the only quantities that are invariant 
when we allow the gap at quasienergy tt/T to close are 
the numbers for fc = 0,7r/a. From similar ar¬ 

guments as above, the expression should be symmetric 
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in fc = 0 and k = tt/ a, and we conclude 
J^ph(O) = 


( 20 ) 


where Qo and Qw/a were the numbers of zone-edge sin¬ 
gularities at crystal momentum 0 and 7r/a, respectively. 
The results in Eqs. (19) and (20) agree with those ob¬ 
tained in Ref. [321 We see that a system can have a non¬ 
trivial edge mode spectrum if the evolution contains an 
odd number of zone-edge singularities, even if the bulk 
Floquet operator is trivial (i.e., if 770 = Vn/a = 0)- 

We finally briefly discuss the special case of two bands. 
For two-band systems, it is always possible to choose 
It/ji) = (1, 0) and \ 1 IJ 2 ) = (0,1) in Eq. (16). Then, the sign 
of the parameter A in Eq. 0 is forbidden to change, and 
defines a conserved charge for each singularity. The net 
charges of zone-edge singularities at A: = 0 and k = tt/ a 
define two topologically invariant Z indices of the system. 
This richer classification for two-band systems relies on 
unbroken translational symmetry, however, and we do 
not expect it to survive disorder. 


C. Time-reversal symmetry with = — 1 


As a final application, we now apply our framework 
to periodically driven systems with time-reversal sym¬ 
metry, identified as a time non-local symmetry of the 
form in Eq. (13) where the symmetry operator S is anti¬ 


unitary and squares to —1 (analogous to AZ class All). 
The condition 5"^ = — 1 implies that the system must 
have an even number of bands, taken to be 2N in the 
discussion below. Following the discussion in Sec. [Ill A[ 
the presence of this symmetry implies that a basis and 
time-origin exist such that the Hamiltonian of the sys¬ 
tem satisfies iJ(k,t) = aH*{—'k,T — t)a, where ct is a 
unitary, Hermitian matrix that is purely imaginary in 
the specified basis (see also, e.g.. Ref. mi). In this basis, 
the evolution operator U{k,t) satisfies 


U{k,t) = aU{-k,T - tyu'^{-k,T)a. (21) 


Note that by substituting t = T, and using 17(k, 0) = 1, 
we recover the time reversal symmetry condition on the 
Floquet operator itself: U{k,T) = aU'^{—k,T)a. 

We begin our discussion of time-reversal invariant 
(TRI) systems below by constructing the general bulk- 
edge correspondence for two-dimensional TRI periodi¬ 
cally driven systems. In doing so, we find that the edge 
mode parities in the bulk gaps can be non-trivial even if 
the system has trivial Floquet bands (see also Ref. m|) . 
This allows driven systems to have edge mode spectra 
that are impossible to obtain in non-driven systems. In 
particular, as we demonstrate subsequently in an illus¬ 
trative example, a driven two-band system can have he¬ 
lical edge modes in its Floquet zone-edge gap. This be¬ 
haviour is in contrast to that of a non-driven two band 
system, where Kramers’ theorem guarantees a gapless 
bulk. As for the cases with no symmetries or with par¬ 
ticle hole symmetry, we find a close connection between 




FIG. 5. Phase bands and topologic al sin gularities for systems 
with time-reversal symmetry, Eq. (211. After deforming to 
a time-periodic evolution U (k, t) which has its zone-edge gap 
open, two distinct classes of evolutions are possible, with an 
even or odd number of pairs of zone-edge topological singu¬ 
larities. a) Phase bands of U{t) for a system with an odd 
number of singularities, and hence one pair of protected he¬ 
lical edge modes its zone-edge gap. b) Phase bands of U{t) 
for a trivial system, yielding no edge modes in its zone-edge 
gap. c) Example distribution of zone-edge singularities, with 
their locations in the t-plane depicted. The charges are in¬ 
dicated by the hlling of the dots, and bars indicate conjugate 
partners, d) Quasienergy band-structure of the model (301, 
when defined in a strip geometry along the j/-direction. Here 
we only display the edge spectrum for one edge. 


these anomalous edge mode phenomena and the appear¬ 
ance of topological singularities in the zone-edge gap of 
the bulk evolution. After working through the example, 
we conclude with the topological classification for three- 
dimensional TRI systems. 


Bulk-edge correspondence for TRI systems in two dimensions 


To simplify the derivation of the bulk-edge correspon¬ 
dence for two-dimensional (2D) TRI systems, we start by 
highlighting some general properties of the Floquet bands 
of such systems. For a periodically-driven system with 
time reversal symmetry in any dimension, Floquet bands 
2m — 1 and 2m are related by time-reversal symmetry 
and the gap between the bands closes at the inversion- 
invariant points in the Brillouin zone (in accordance with 
Kramers’ theorem). This holds for all m, and only “even” 
gaps 2m may thus remain open. Here we use the specific 
assignment of even and odd indices defined by the order¬ 


ing scheme of Sec. IB 


In a finite geometry, a 2D TRI system can have pro¬ 
tected edge modes in its bulk gaps. Time reversal sym- 
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metry requires any chiral modes to come in time-reversal 
conjugate pairs; Kramers’ theorem guarantees that an 
odd number of such pairs cannot be gapped out by any 
time reversal symmetry preserving perturbation. The 
parity i^tr(w) of the number of such “helical” edge mode 
pairs appearing in gap 2m is thus a topological invariant 
that can be associated with gap 2m. 

We now set out to find an expression for j^tr(w)- 
First, noting that for each m the Floquet bands 2m and 
2 m — 1 are related by time-reversal symmetry, to each 
such pair we associate a Fu-Kane Z 2 index Zm, just as 
for the bands of a TRI non-driven systenP^l. Spectral 
flow argument^^, which must hold for both static and 
Floquet bands, show that the relative edge mode parities 
of gaps 2n and 2n — 2 (i.e., above and below the pair of 
time-reversal conjugate bands 2n and 2n — 1) are cap¬ 
tured by Zn- For a system with 2N bands, we can thus 
write 

m 

I^Tr(w) = tZTR{N) Zn, (22) 

n—1 


where J^TR(fV) is the edge mode parity in the zone-edge 
gap. Note that nr!,=i umst always be unity; the pref¬ 
actor ^tr(-N) ensures that the correct edge mode parity 
is recovered when setting m = N in Eq. ( [2^ . 

To find we first simplify the symmetry condi¬ 

tion in Eq. (21) by smoothly deforming C7(k, t) into an 


evolution whose Floquet operator is the identity. During 
the deformation we preserve TRI and keep the zone-edge 
gap open, such that the edge mode parity vtr{N) in the 
zone-edge gap remains unchanged. An example of such 
a deformation is given in Ref. HU We first define an ef¬ 
fective Hamiltonian Ffeff(k) = P„(k, r)e„(k), where 

P„(k, T) is the Floquet eigenstate of band n, and e„(k) = 
(/)„(k, T) /T is the corresponding quasi-energy of the band 
(see Sec. |IB| for the quasienergy zone convention). The 
family of evolutions C/(k, t;a) = [/(k, for 

a € [0,1] defines a smooth interpolation of evolutions 
from U to a time-periodic evolution U, given by 




(23) 


The Floquet operator associated with the evolution 
C/(k,t;a) has quasi-energies {(1 — a)en(k.,T)/T}-, the 
spectrum uniformly contracts, and thus the zone-edge 
gap stays open throughout the interpolatiorl^. Given 
that TRI is also preserved during the deformation, 
i^tr(A^) must be the same for U and U. 

Using the fact that C/(k, T) = 1, symmetry condi¬ 
tion (21) simplifies to U(k,<) = aU*{—'k,T — t)a. An 
evolution operator with this property can be written as 


N 

U{k,t) = ^ [p„(k,t)e-*^"(‘^’*) +P^(k,t)e-*^-('^’‘)] , 

71 — 1 

(24) 

where the bands {P„, Pfi} for n = 1... A^ are all orthogo¬ 
nal. We give a specific prescription for defining the labels 
n and n below. 


In contrast to the particle-hole symmetric case, TR 
symmetry relates conjugate bands at different times: 

Pnik,t)=aP^i-k,T-t)a (25) 

(j)nik,t) =-(j)n{-k,T-t). (26) 

For the last equality we used that ^„(k, T) = 0. We 
assign labels to the bands such that all (j)n{k,T/2) are 
positive, and then order these bands according to in¬ 
creasing (fin- Their conjugate partners then follow from 
Eq. (26). In Figs. ,b we show two examples of time- 


periodic phase-band structures that satisfy the time re¬ 
versal symmetry above. 

We now smoothly deform one half of the phase bands of 
U, via the phases {(/)„(k, t)}, in a way that preserves the 
boundary condition 0„(k, T) = 0, while keeping the pro¬ 
jectors {P„(k,f)} constant. The other half of the bands 
follow in accordance with the symmetry above. Through 
considerations similar to those made in Sec. we find 
that the phase bands of U can be deformed to zero every¬ 
where in k, t-space except around its zone-edge topolog¬ 
ical singularities (which connect bands N and N). The 
edge mode parity i^tr(A^) is thus completely determined 
by the constellation of zone-edge singularities of P. 

Next we consider which features of the zone-edge singu¬ 
larity constellation are topologically invariant. As for the 
case without symmetries, through smooth deformations 
of U we can annihilate singularities of opposite charges. 
If U has a zone-edge singularity with charge q at (k, t ), 
symmetry dictates that it has another with charge —q 
at (—k, T — t); see Figj^ for an example distribution of 
zone-edge singularities. 

In general, when two singularities annihilate, their con¬ 
jugate partners must annihilate as well. Importantly, 
conjugate singularities cannot annihilate directly with 
each other (see Appendix]^. Therefore it is only possi¬ 
ble to annihilate two singularity pairs at a time. To give 
a concrete example, if singularities 1 and 2 in Fig. an¬ 
nihilate each other, then by symmetry singularities 1 and 
2 will annihilate as well. Singularities 3 and 3 will then 
remain, with no way to be eliminated without closing the 
zone edge quasienergy gap. 

According to the arguments above, the parity p of the 
number of singularity pairs is invariant under any smooth 
deformation of the evolution that keeps the zone-edge 
gap open. Conversely, any two evolutions U, U' with 
the same parity p can be deformed into each other. In 
particular, if U has an even number of singularity pairs, it 
can be smoothly deformed into the identity. Evolutions U 
therefore fall into two classes: those with an odd number 
of zone edge singularity pairs (Fig. [^) and those with 
an even number of pairs (Fig. §>). Evolutions within the 
same class must have the same edge mode parity !ztr(A^), 
since they can all be related by smooth deformations that 
preserve ^tr(A^)- Thus zztr(A^) must depend only on 
(—1)^. Since we already found that evolutions with even 
p are topologically trivial (i.e., smoothly deformable to 
the identity), we identify 

J^tr(A^) = (—1)^- 


(27) 
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The edge mode parity in the zone-edge gap is thus given 
by the parity of the number of singularity pairs of U. 

Combining results (22) and (27), we then obtain for 
the edge mode parity of gap m: 


VTK{m) = (-1)^ Zn, (28) 

n—1 

where p is the parity of the number of zone-edge sin¬ 
gularity pairs of C/(k, t). We see that even if the bulk 
Floquet operator of a system is trivial (i.e., = 1 for 

all n), the edge mode spectrum can be nontrivial. This 
is precisely the case when 17(k, t) has an odd number of 
singularity-pairs in its zone-edge gap (see Fig. [^b). 

We end our discussion by noting that the “anoma¬ 
lous” Z 2 index = (—1)^, Eq. ( [^ , can be 

found directly as a time-domain expression in terms of 
C/(k, t). According to the time-reversal symmetry condi¬ 
tion t/(k, t) = aU*{—'k,T — t)a, one singularity out of 
each conjugate pair must occur within the first half of 
the driving, 0 < t < T/2. Therefore p is in fact equal to 
the number of zone-edge singularities of V {t), where V (t) 
is the evolution given by U (t) restricted to the first half 
of the drivin^^. In Sec. |ll] we found that the number of 
zone-edge singularities of an evolution is given by wind¬ 
ing number of the evolution in its zone-edge gap. Hence 
we have 


P = W 2 iv[E(<)] (mod 2). (29) 

The above result is consistent with the results obtained 
in Ref. l44l 


Example: Nontrivial Z 2 index for a two-band TRI system 


To illustrate one of the new topological phenomena 
which arise in periodically driven systems, we now show 
that periodic driving allows a 2D TRI system with two 
bands to have protected helical edge modes. This is in 
contrast to the situation for non-driven systems, where a 
minimum of four bands is required. We explicitly demon¬ 
strate this behavior for a specific model, using the results 
above as well as direct numerical calculation. 

Consider a spin-1/2 particle on a square lattice with 
one orbital per site. We construct the evolution based on 
the model with non-trivial winding numbers presented in 
Sec. m now with spin playing the role of what was an 
orbital index. Specifically, we define a TRI evolution via 
the time-periodic Bloch Hamiltonian 


H{k,t) 


Ho{k,2t), t<T/2 

ayHS{K2{T-t))ay, t>T/2 


where, as a 2 x 2 matrix, iJo(k, t) has the same form 
as the Hamiltonian presented in Sec. |H B| The Hamilto¬ 
nian H{k.,t) is by construction time-reversal symmetric. 
In real-space, Hamiltonian (30) contains on-site terms 


as well as nearest-neighbour and next-nearest neighbour 
hopping terms. 

In Ref. 03] it was noted that the Floquet operator 
t/o(k, T) associated with i7o(k, t) is unity for the param¬ 
eter values 5ab = 0, A„ = 57r/(2T), and that the cor¬ 
responding time evolution operator Uq (k, t) has winding 
number 1. For our model ( |30| ), we choose the parameters 
of Hq such that 6ab = 0 and A„ = fi-n jT. The Flo¬ 
quet operator of the (translationally invariant) system 
governed by H is again equal to the identity. 

We now use Eq. ( [^ to predict the number of helical 
edge mode pairs, = 1), that will appear for this 

model when defined in a strip geometry. First, note that 
for the system without edges the Bloch-space evolution is 
periodic with its zone-edge gap open, and hence U (k, t) = 
C/(k, t). Second, because the evolution in the first half of 
the dri ving corresponds to that of the model discussed 
in Sec. HB we know that the evolution for 0 < t < r/2 
possesses exactly one zone edge singularity. Hence p = \ 
and we have 


^^tr(1) — (—1)^ — —1- (31) 

We thus expect one helical edge mode pair to appear in 
the zone edge gap, for each edge. 

To confirm the analysis above, we numerically study 
this model in a strip geometry, with edges parallel to 
the y-direction. The Hamiltonian of the strip geometry 
system is obtained by truncating the real-space tight- 
binding Hamiltonian in the ^-direction. Using this trun¬ 
cated tight-binding Hamiltonian we numerically calcu¬ 
late the corresponding Floquet operator, and obtain the 
quasienergy band structure as a function of the conserved 
momentum component ky, see Fig. At each edge we 
find a pair of helical modes, which are time-reversal con¬ 
jugates of each other. This model thus explicitly demon¬ 
strates the existence of anomalous helical edge modes in 
two-dimensional periodically driven systems with TRI, 
and shows that non-trivial topology can be found even 
in a case with only two bands. 


time-reversal symmetry 


The Z 2 index for two-dimensional systems, see 
Eqs. (28) and (29), can also be used to define a Z 2 index 
i^ 3 d(w) for three-dimensional (3D) systems with time- 
reversal symmetry. The index iy 3 u(jn) indicates whether 
or not Floquet gap m hosts non-trivial surface states, 
analogous to those of a strong topological insulator. In 
the same way as for the non-driven case, we consider the 
two 2D TRI systems defined by the restriction of U (k, t) 
to the planes ki = 0 and ki = tt / o , where i can be x,y 
or z. We then identify the index as the product 

of the indices i^tr(r^), from Eq. (28), for the two sys¬ 
tems. Note that if the indices calculated for the planes 
ki = 0 and ki = tt/o are the same, their common index 
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determines whether the system acts as a weak topological 
insulator in the plane orthogonal to ki. 


IV. DISCUSSION 

In this paper we found that the “phase-band struc¬ 
tures” of time evolution operators provide a powerful 
basis for visualizing and understanding the topology of 
Floquet-Bloch systems. By considering smooth deforma¬ 
tions of the phase bands, we showed that topologically 
protected degeneracies, or topological singularities, play 
a crucial role in distinguishing the topology of driven and 
non-driven systems. In particular, the presence of phase 
Brillouin zone edge topological singularities can present 
an obstruction to smoothly deforming the evolution of a 
driven system into one obtainable in a non-driven system. 

Our work demonstrates a general method for topo¬ 
logical classification of Floquet-Bloch systems, based 
on identifying all features of the phase band structure 
of a given system that cannot be removed by smooth 
deformations. This approach appears to offer means 
to exhaustively classify Floquet-Bloch systems and to 
straightforwardly derive the corresponding bulk-edge cor¬ 
respondences. 

In the cases we considered, with symmetries imposed 
in analogy with the ten Altland-Zirnbauer classes, we 
found that the edge spectra associated with individual 
bulk gaps of the Floquet operator have the same fea¬ 
tures as those of non-driven systems in the correspond¬ 
ing symmetry classes. However, we found in all cases 
that periodic driving could induce global edge spectra 
that are impossible to obtain in non-driven systems. In 
particular, with periodic driving, topologically protected 
helical edge states can be produced in time-reversal in¬ 
variant systems with only two-bands, while a minimum 
of four bands is needed without driving. In each case we 
considered, these new “anomalous” features were closely 
related with the appearance of zone-edge singularities in 
the time-bulk. These phenomena further demonstrate 
that the relation between the topological properties of 
the bulk evolution and the appearance of protected edge 
modes is fundamentally changed in the driven context: 
the topology of a periodically driven system cannot be 
fully characterized by the stroboscopic Floquet operator 
U (k, T) or the corresponding effective Hamiltonian alone. 

Finally, our results provide new intuition about the 
topology of Floquet-Bloch systems. While time-domain 
invariants such as the winding number found in Ref. |3SI 
offer a mathematically well-defined prescription for char¬ 
acterizing topology in driven systems, often a clear phys¬ 
ical picture is missing. Here, one of our central re¬ 
sults is that non-trivial topological phenomena in peri¬ 
odically driven systems appear when topological singu¬ 
larities are introduced into the phase bands of the bulk 
time-evolution operator. In particular, in any driven sys¬ 
tem where the Floquet bands have different topology 
from that of the initial instantaneous Hamiltonian (i.e., 


at t = 0), at least one topological singularity must be 
encountered during the evolution. We expect that this 
insight may help provide guidance for the construction 
of new driving protocols to realize non-trivial topological 
phases in periodically driven systems. 

The instantaneous and time-nonlocal symmetries con¬ 
sidered in Sec. |HIA| were chosen to illustrate the power of 
the phase band framework. The symmetries altered the 
topology of the bulk evolution through their ability to 
protect new types of topological singularities. However, 
we expect that other types of conditions could give rise to 
new non-trivial topology in the phase bands. The explo¬ 
ration of other symmetries and their role in the topology 
of periodically driven systems, as well as the search for 
good candidate systems for their realization, are interest¬ 
ing directions of future study. 

We gratefully acknowledge E. Berg and N. Lindner for 
enlightening discussions. MR thanks the Villum Founda¬ 
tion and the People Programme (Marie Curie Actions) 
of the European Unions Seventh Framework Programme 
(FP7/2007-2013) under REA grant agreement PHF-GA- 
2013-627838 for support. 


Appendix A: Time-domain expressions for the 
invariants vg and Vi 

In this appendix we derive a time domain integral ex¬ 
pression for the bulk invariant WmiU] in Eq. Q, which 
corresponds to the number of edge modes appearing in 
gap m when the system is defined in a geometry with 
edges. We directly show that this invariant is equivalent 
to the winding number invariant IU[C4] of Ref. [43l with 
e set equal to a quasienergy value inside gap m: 

/ rl'^krlf 

^Tr {UldtUm [Uldk^U^, Uldk,U ^]} . 

(Al) 

Here Um is a time-periodic evolution operator, satisfy¬ 
ing Um{^,T) = 1, which is obtained from the original 
evolution U by a smooth deformation in which gap m of 
the quasienergy spectrum is kept open. Explicitly, the 
time-periodic evolution operator Um can be obtained by 
deforming the final values of the phase bands 0„(k, T) to 
zero for n = 1.. .m, while the final values of the remain¬ 
ing bands are deformed to 27r. 

In order to demonstrate the equivalence of the two in¬ 
variants, i.e., to show Wm\U] = IU[C4n], we first consider 
two special cases and then discuss the general situation. 
In this appendix we label points in the three dimensional 
k, t-space by a single dimensionless vector 

k t \ 

27r/a’ ry ■ 
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at the edge of the neighbourhood to tt and — tt at the cen¬ 
ter (letting (j)d go continuously to tt in the process), see 
Fig.i Under the deformation we keep the eigenstates of 
Umis) constant everywhere. 

The deformed evolution operator Um is equal to the 
identity for all s, except in the small region of radius Ss 
that surrounds the singularity (see Fig. |^. Within this 
neighbourhood, Um takes the form 


FIG. 6. The deformation of an isolated singularity discussed 
in Sec. A 2 In the figure, phase band 1 (blue) is shifted up 
by 2n for clarity of illustration. 


1. Winding number in case of no singularities 

To begin, we first consider the case where Um (k, t) has 
no topological singularities in the zone-edge gap. It is 
then possible to continuously deform the evolution op¬ 
erator Um into one corresponding to a non-driven sys¬ 
tem, as described in Sec. [l] In doing so, the winding 
number W[Um], a topological invariant, cannot change 
its value. After the deformation, W[Um] is simply the 
winding number of a non-driven system described by the 
Hamiltonian 


i^m(k) = — ^ P„(k,T). (A2) 

The winding number of a system governed by such a 
Hamiltonian was found in Ref. |43]to be Cn, where 

Cn is the Chern number of Floquet band n. Hence 

m 

IF \Um\ = ^ Cn, (no zone-edge singularities). (A3) 

n—1 

Note that the winding number will always be zero in 
the quasienergy zone-edge gap if the phase bands do not 
host any zone-edge singularities. This follows from the 
fact that the sum of Chern numbers for all bands must 
evaluate to zero, — 0- 


2. Winding number in the case of one singularity 

We now consider the case where (/>„(k, T) = 0 for all 
n, and Um (k, t) has only one singularity in the zone-edge 
gap, located at Sg = (ko/(27r/a),tg/F). At the singu¬ 
larity, the two touching bands N and 1 have phases 4>d 
and (j)d — Stt, where (pd is a real number determined by 
details of the evolution. In this case we can deform the 
phase bands to zero for all k, t, except in a small spherical 
neighbourhood of radius Ss that surrounds the singular¬ 
ity (here lengths are computed with respect to the usual 
norm on the dimensionless vector s). Within the neigh¬ 
borhood, all Af — 2 bands not involved in the singularity 
can still be flattened. The phase values of the two inter¬ 
secting bands are deformed to evolve linearly from zero 


t^m(s) = \Xn){Xn\ + ^ | V'a) (V'b I (s), (A4) 


n^l,N 


a,6=1.AT 


where Mab{s) is a 2 x 2 matrix whose eigenvectors are 
the eigenvectors of the matrix (s — SQ)jSjk(Jk- Here, the 
real invertible 3x3 matrix S was defined in Sec. m 
and {cTfe} are the Pauli matrices. From the description 
of the flattened phase bands above, we know that the 
logarithms of the eigenvalues of M must grow linearly 
from 0 at |s — Sg| = (5s to —ztt and Itt at s = Sg. For 
|s — Sol < 6s, Mab thus takes the form 


Mnbis) 

+ i;+(s)i;+*(s)e-*’^(l"-"«l/<5«-i). 


(A5) 


The vectors v~(s) and i;“''(s) are the eigenvectors of the 
traceless 2x2 Hermitian matrix (s — Sq) jSjkcrk, cor¬ 
responding to negative and positive eigenvalue, respec¬ 
tively. Since the matrix S is real and invertible, we can 
write it as S' = R 1 AR 2 , where i?i and R 2 are orthogonal 
and A is a diagonal matrix with positive entries (this is 
the singular value decomposition of S, see e.g.. Ref. I351l . 
A continuous deformation of the entries of A to 1 re¬ 
sults in a orthogonality-preserving continuous interpo¬ 
lation of the eigenvectors ^^(s) to the eigenvectors of 
Rjki^ — So)j(Tfc, where R = By continuously de¬ 

forming the vectors in this way, Um is deformed into 
an evolution operator Vm still of the form (A4), but with 
the matrix M given by 


M(s) = 


1 , 


|s — Sg| > 6s 


-exp [^(s - Sg),%(Tj-] , |s-So| < 5s. 

(A6) 

Recall that R is orthogonal and its determinant |i?| 
is the charge of the singularity, q = |i?i||i? 2 | = |S|, see 
Eq. (|^. In Appendix [b] we explicitly evaluate the wind¬ 
ing number of the evolution operator Vm- We find: 


W[Vm\ = -\R\ = -q. 


(A7) 


Using the fact that the winding number could not change 
during the deformation from Um to Vm, we thus establish 


W[Um] = -q. 


(AS) 


In other words, if Um contains one isolated singularity, 
the winding number of Um is given by the corresponding 
charge of the singularity (with a minus sign). 
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3. The general case 


We now consider the general case, where Um has 
N topological singularities in the zone-edge gap, with 
charges In order to evaluate the winding number, 

we deform [/^(k, <) as described in Sec. IIA and shown 
in Fig. The deformed evolution is the identity every¬ 
where except for in small isolated regions surrounding the 
singularities (region I), as well as in the short ramping 
time-interval St at the end of the driving (region II). 

The winding number ( |Al[ ) is defined as an in¬ 
tegral over k, t-space of the quantity Fm (k, t) = 
^FT{Ul,dtUrn[U^rndk^Um,Ul^dk^Urn])- For the de¬ 
formed, “band-flattened,” system, is only nonzero 
in each of the isolated regions that surround the singu¬ 
larities, and in the final ramp region II. We can therefore 
split up the integral of Fm into a sum of integrals over 
each of these nontrivial regions. 

From the first special case we examined, i.e., for an evo¬ 
lution with no singularities, we know that the integral of 
Fm(k,t) over region II equals where {C„} are 

the Chern numbers of Floquet bands 1... to. From the 
second special case, we know that the integral of Fki(k, t) 
over one of the regions surrounding a zone-edge singular¬ 
ity equals minus the charge of the singularity, i.e., —q. 

Summing the integrals over all regions, we obtain: 


N 


W[Um] = Y.C^-Y. 


Fi 


(ra) 


(A9) 


n—1 


where Sijk is the Levi-Civita symbol. Summation over 
repeated indices is used and will be used in the rest of 
this appendix. 

In order to exploit the s-space spherical symmetry of 
the deformed evolution Vm, we shift from Cartesian co¬ 
ordinates to spherical coordinates centred around Sq, de¬ 
fined such that 


s — So = (s sin0sin(/), s sin0cos0, s COS0). (B2) 


After the coordinate transformation, W is expressed as 

i-Ss p-n p2Tr 

IF = ^ y dsjdej dct>\J\J,^JjpJkj 

■ Tr {FPd^M M'fdpM M^d^M} , (B3) 

where J is the Jacobian matrix of the coordinate trans¬ 
formation, and the Greek letters run over the co¬ 

ordinates s, 0, (j). We now use the following useful identity 
for the Levi-Civita symbol that holds for any real invert¬ 
ible 3x3 matrix 


-^ia. 0-^k'y^ijk 




(B4) 


With the help of this identity we see that the Jacobian 
matrices always cancel out: 


W = 


^a /37 

247r2 




ds / dO 


d(f) 


TrlM^d^MM^dpMM^d-yM} . 


Note that Um can be constructed by deforming the phase 
bands of U only at the end of the driving. Therefore the 
net charge of all zone-edge singularities in the time-bulk 
should be the same for U and Um- Thus = 

where are the zone-edge singularity 

charges for the original system with evolution governed 
by U. Hence, we finally have the result for the number 
of edge modes in a two-dimensional system: 

m 

W[Um] = E - E 

n—1 i 

This is what we set out to show. 


Summing over the indices, we obtain: 


W = 


1 

Stt^ 


ds / dd 


d(j) 


Tr {M‘^dsM [MiSeM, . 


Using the cyclic property of the trace as well as the iden¬ 
tity dMM^ = —MdM\ we get 




w = - 




ds / dO 


d(f) 


Tr {M'^dsM [d 0 M\ d^M]} . (B5) 

We now consider the explicit canonical form of M, 
Eq. (A 6 ), in polar coordinates, in the region |s —So| < 6 s: 


Appendix B: Derivation of Eq. (A 81 


In this appendix we prove that the winding number 
of an evolution operator Vm of the form in Eq. ( |A4[ ), 
with the matrix M given in (A 6 ), is equal to —\R\- We 
begin by inserting Vm from Eqs. (A4| and (A 6 ) into 
Eq. (lAll), to obtain 


ITT , 


M(s, 6 »,((i) = -exp |^- — (s-So)-rj , (B 6 ) 

where = RijCrj. Defining s{9,4>) = (s — So)/s, we 
evaluate each of the factors in the integrand 


M^dgM = — ^ s ■ T 
ds 


f d^sTr{M^ds,MM^ds^MM^ds^M} , 

J \s—sq\<Ss 

dsM = zsin 1 

nTs\ 

(Bl) 

d^M = i sin 1 

^Js) 
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Hence, after performing the integral over s, we obtain 

W = [ d,9 [ d(/)Tr{s • r [Sgs • t, d^s ■ r]} . (B7) 

IGtt Jo Jo 

Working on the integrand, we note: 

Tr{(s-T)(aes-T)(a,^s-T)} 

= Sj deSj dtpSk RiaRjbRkc Tr {(TaabCTc} ■ 

(B8) 

Using the Pauli matrix identity Tr{(TiCTjtTfe} = 2ieijk, we 
obtain 


Tr{(s-T)(aes-r)(50S-r)} = s^dgSjd^Sk ■ 2iR iaRjbRkc^abc 

* 2z|-/?|£:2j7c' 

Here we used the Levi-Civita symbol identity ( |B4[ ), and 
the fact that = |i?|. Restoring antisymmetry in 0 

and (j), and going back to vector notation, we have 


Tr{(s-r)[(aes-T), (a^s-r)]} = B|i?| s ■ {dgs x d^s) 

= B|i?|sin0. (B9) 


Hence the integrand in Eq. (B7) is simply 4i|i?| times the 
surface area element of the sphere. We thus have 


where both the projectors Pn{t) and the phases are 
continuous functions of t, and we suppress the momen¬ 
tum index for brevity. Here, the projector is the 
transpose of P„ in the basis where H{k, t) = —H*{—k, t). 

When we impose the requirement that the projectors 
{Pn{t)} are continuous in t, the labeling scheme defined 
in Sec. |HI B| cannot be employed. In particular, in the 
original scheme the phases {(j)nik,t)} were constrained to 
live within the interval [0,7r], “bouncing” off the phase 
Brillouin zone edge when meeting a conjugate partner 
band there. Under the new construction, continuity of 
the projectors requires that the phases of conjugate bands 
pass through each other, such that a phase (pnit) may 
move from inside to outside the interval [0,7r]. We note 
that such crossings of partner bands are protected by 
symmetry and can occur at phases which are either even 
or odd multiples of tt. We identify them respectively as 
phase-0 and phase-7r topological singularities. 

We now use this picture of topological singularities to 
show that the Pfafhan of i7eff changes its sign when a 
phase-TT singularity passes through the plane t = T. We 
first express the effective Hamiltonian at the inversion in¬ 
variant point fcinv, Heff(A:inv), in terms of the continuously 
defined phase bands (/>„ and 


d9 / d(^ 4i|i?| sin0 

Jo 

= -\R\ 

Recalling that q = sgalS*] = |i?|, we see that the winding 
number contribution of an isolated singularity is given by 
the charge q of the singularity. 



Appendix C: Sign change of the Pfafhan for 
zone-edge singularities passing through t = T 

In Sec. im of the main text, on the topological clas¬ 
sification of one-dimensional systems with particle hole 
symmetry, we made use of the fact that the sign of 
Pf [iJeff] switches when a phase- tt (i.e., zone edge) singu¬ 
larity passes through the plane t = T due to a continuous 
deformation of the evolution. In this appendix we give a 
proof of this statement. 

We begin by introducing an alternative labeling for 
the phase bands at the inversion-invariant points in the 
Brillouin zone. A crucial point in this paper is that it 
is not always possible to define the eigenstate projectors 
{Pn{k,t)} of the evolution operator such that they are 
continuous everywhere in (fc, t)-space. However, we can 
always define them such that they are continuous along 
an arbitrary line in /c,t-space. At an inversion-invariant 
point fcinv in the momentum Brillouin zone, we can there¬ 
fore write the evolution operator of a particle-hole sym¬ 
metric system in the form 


n 


(CI) 


n 


Pu{T) 



(C2) 


Here F is the modulo 27r function, defined such that it 
takes value between —tt and tt. The Pfafhan of iJeff can 
then be written as (see, e.g., Ref. [55|: 


Pf [i7eff(fcinv)] = Pf 


E(^n(r)-P„^(T)) 


n 

(C3) 

Above we identihed phase- tt singularities as crossings of 
partner bands (l)n(t), through the phase Brillouin 

zone edge at phase tt. Hence, under a continuous de¬ 
formation where a phase-7r singularity passes through 
the plane t = T, the value of 4>n{T) for some n crosses 
through an odd multiple of tt. As j)n{T) passes through 
the quasienergy zone edge, P(0„(T)) therefore changes 
it sign. During this process, F{(j)n{T)) for all other n and 
all of the eigenstate projectors {Pn} may only change in- 
hnitesimally. From the above expression for the Pfafhan, 
Eq. (C3), it then follows that Pf [i7eff(fcinv)] changes its 
sign when a phase- tt singularity passes through t = T. 


Appendix D: Proof that TR-conjugate singularities 
cannot annihilate 


In this appendix we show that two TR-conjugate zone- 
edge singularities cannot annihilate with each other, as 
claimed in Sec. EH ^ The argument follows indirectly 
from the main text. 

The discussion in Sec. |HI Cj implies that if it were possi¬ 
ble for conjugate singularities to annihilate, the evolution 
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of any system with TRI could be smoothly deformed to 
the identity, while keeping the zone-edge gap open. How¬ 
ever, in Sec. |III C| we saw an example of a system with one 
pair of helical edge modes in its zone-edge gap. Suppose 
it were possible to deform the evolution of this system 


into the identity, without breaking TRI or closing the 
zone-edge gap. Then the helical edge modes would have 
to disappear without a quasienergy gap closing. However, 
this is impossible and implies a contradiction. Hence it 
is not be possible for a pair of time-reversal conjugate 
singularities to annihilate each other. 
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